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We consider the time, angular and energy distributions of SN1987A events discussing the quality 
of their agreement with the expectations. A global interpretation is performed considering a simple 
model, based on the standard scenario for the explosion. Despite the contrasting and confusing 
indications, a straightforward fit to the data provides a result that does not contradict but rather 
supports the expectations. The calculated electron antineutrino flux is applied to predict the relic 
neutrino signal. 



I. INTRODUCTION 

The Symposium held in Moscow on February 2007 |l| 
was devoted to discuss open problems of SN1987A. The 
meaning of LSD events [2] was discussed, and the po- 
tential of the model of Imshennik and Ryazhskaya Q 
was shown. The observations of Kamiokande-II |4j , 1MB 
Q and Baksan @ were covered by the contribution of 
Alexeev. In the spirit of the Symposium, we recall that 
(1) certain features of these observations are difficult to 
understand, see e.g. 0j; (2) there is not yet a complete 
model for the explosion, that could shed doubts on the 
assumption that we know how neutrinos are emitted; this 
applies also to ref. [3| . 

In this work that expands the talk of F.V. at [1[ we 
study the features of Kamiokande-II, 1MB and Baksan 
data and analyze them in the 'standard scenario' of emis- 
sion, first advocated by Nadyozhin 0]-also called Bethe 
and Wilson scenario delayed explosion or ^-assisted 
explosion. We recall its features: On top of a several- 
seconds lasting emission from the proto neutron star, the 
cooling, there is an initial intense emission lasting a frac- 
tion of a second. The v e and v e irradiated in this phase 
transfer energy and help the stalled shock wave to re- 
sume, eventually leading to explosion. We call this phase 
of neutrino emission accretion. For progresses toward 
the implementation of the 'standard scenario' in numer- 
ical simulations, see [Icj . 

II. FEATURES OF THE DATA 

The observed events are 11 or 12 in Kamiokande-II, 
the 6 th being below the threshold of 7.5 MeV. The du- 
ration and the energy of the signal is not safely known. 
Thus we consider a unified time window of T = 30 s in 
this work. The candidate signal events increases to 16 
in Kamiokande-II (the events number 13,14,15,16 being 
below the threshold), it is 8 in 1MB and 5 in Baksan. Re- 
call that most signal is due to v e p — > e + n [IBD reaction] 
that has the largest cross section; unless stated otherwise 
we assume that all signal events are due to IBD. 

The background in Kamiokande-II declines rapidly 
with the number of hit phototubes; on average, there 



should be about 0.7 background events above threshold 
and 5.6 in total. Kamiokande-II analyzed the whole vol- 
ume of the detector. The region close to the walls (in 
particular the upper wall) is the less safe against the risk 
of background. Five of the 11 events are close to the 
walls, lying in the outermost 4% of the volume of the 
detector; they include the events number 3,4,10 of the 
dataset, close to the threshold. We expect 1 background 
event in Baksan in T = 30 s; this, together with a typi- 
cal expected signal of about 1.6 antineutrinos (ju = 2.6) 
gives a reasonable Poisson chance /i™e~ M /n! = 8 % when 
n = 5. 1MB instead can be assumed to be background 
free. 

Energy distribution This distribution is difficult to in- 
terpret. We select 3 questions for the discussion: 

(i) Do the observed energies meet the expectations? 

(ii) Are the average energies of 1MB 31.9 ± 2.6 MeV and 
of Kamiokande-II 15.4 ± 2.4 MeV compatible? 

(iii) Why the first 4 Kamiokande-II events have average 
energy 12.6±2.4 MeV and the last seven 17.1 ±3.3 MeV? 
The first can be answered fitting the data with a black 
body (thermal) v e spectrum with luminosity oc R^.T^. 
We confirm that the Kamiokande-II and 1MB datasets 
are compatible only when we consider the 90% C.L.; Bak- 
san is compatible with both. When we combine the data 
in a unique analysis, the values of the radius is larger 
than expected, R c ~ 30 km. The second question can 
be discussed quantitatively calculating the average ener- 
gies of the events, keeping the antineutrino temperature 
T c as a free parameter (see e.g., [HI): Kamiokande-II 
and 1MB suggest different values of T c . The first two 
questions are the traditional formulations of the "energy 
problem". They suggest two possible ways-out: a) The 
energy distribution is st rong ly non-thermal. This is ex- 
plored, e.g., in [T3| and [l3[. Our objection is that this 
is a major deviation from the expectations that should 
not be admitted in a conservative analysis, as, e.g., [l4T |. 
b) Another possibility is that there are some low energy 
background events in Kamiokande-II [l5| , see the feature 
at E = 4 - 8 MeV in Fig. [TJ but recall that the back- 
ground events above the nominal threshold are expected 
to be few. The third formulation of the problem [l7|],0 
offers a new clue: the excess of low energy Kamiokande- 
II events is due to the early detected events. This is fine 
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FIG. 1: Dots: observed energies in Kamiokande-II in 30 s. Con- 
tinuous curve: expected cumulative spectrum (background + signal 
as in Eq. \3jj as a function of the energy in Me V. The Smirnov- 
Cramer-Von Mises (SCVM) 0/ goodness of fit is 19%; compare 
also with \H\j . 



since the standard scenario predicts the existence of the 
initial phase of accretion possibly with peculiar features. 
The high luminosity helps account the early events seen 
by Kamiokande-II. 

Angular distribution The data and the expected an- 
gular distribution have been recently compared in 
keeping into account the angular bias of 1MB [5]. The 
SCVM test shows that the "problem" of the forward 
peak is not severe: if all event are IBD, the good- 
ness of fit is 6.4% for 1MB and 8.6% in Kamiokande-II. 
With 1 or 2 elastic scattering events (0.3-0.6 expected in 
Kamiokande-II) this improves further. Note that, even 
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FIG. 2: Expected 2- dimensional (positron- energy in MeV, cosine 
with the direction of SN 1987 A) -distribution for Kamiokande-II. 
The contour plot with equally spaced heights shows the mild direc- 
tionality of IBD reaction and the peak of ES events (forward region 
near 15 MeV) whose shape is dictated by instrumental effects (see 
e.g., E2/j. 
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FIG. 3: In the abscissa, the time of the event; in the ordinate, its 
progressive number. Two time distributions, comprising 22 signal 
and 7 background events, are given. Also indicated their goodness 
of fit values. In the lower one, all signal events belong to cooling, 
where the signal declines with a time constant r c = 4 s. In the 
upper one, only 13 signal events belong to cooling; the remaining 
9 belong to accretion phase and are distributed in the interval t < 
0.7 s. The dotted vertical line marks the end of the first second. 



being ready to consider something exotic, it is hard to 
imagine a reaction that is forward-peaked but too much, 
as needed to locate half of the 1MB events in the region 
30° < 8 < 60°. Thus the discrepancy with the expec- 
tations is not very compelling. We also estimated the 
angular distribution taking into account elastic scatter- 
ing events (ES). Assuming oscillations with normal mass 
hierarchy and increasing non-electronic neutrinos as al- 
lowed in [l8| we calculated the time-averaged distribu- 
tion of Fig. [2 From this, we find that the probability 
that the first event detected by Kamiokande-II is due to 
ES is about 30%. 

Time distribution The time sequence of the 
Kamiokande-II events is at first sight odd; there is 
a cluster of six events in the first second and a "gap" 
of 7 seconds between the event number 9 and the 
event number 10, difficult to understand on physical 
basis. One way-out is to combine all data in a single 
dataset as in Fig. [3J This synthetic dataset can be 
compared with the expectations. Besides the known 
background component, n\> = 6.6 expected events, the 
signal should accumulate in two phases: a first one with 
n a events, followed by a second and longer phase with 
n c = 29 — rc;, — n a events: 



N(t) = n b t+n a [l+6{t-T a ){± - l)]+n c [l- 



-t/T e 



(1) 



where the suffixes b, a, c mean background, accretion, 
cooling. The goodness-of-fit values shown in Fig. [3] sug- 
gest that a two component model of this type has no 
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anti-v e energy [MeV] 

FIG. 4: Continuous curve: instantaneous u e flux for M a = 
0.15 Mq and T a = 2.5 MeV. Dotted curve: black body distribution 
with the same luminosity (1.1 X 10 53 erg/s) and average energy 
(IS MeV), i.e., with parameters R c = 82 km and T c = 4.1 MeV. 

problem to reproduce the data. 1 See [13, HH for further 
useful discussion of the time sequence of the events. 



III. A TWO COMPONENT MODEL FOR v„ 
EMISSION 

Let us summarize the previous section. There is some 
disturbing feature of the data. The presence of few elastic 
scattering events does not seem to be essential to explain 
the angular distribution, that conforts us to proceed with 
the usual IBD hypothesis for the signal. We showed that 
background events should be accounted for and that it is 
important to consider very seriously the existence of an 
initial phase of high neutrino luminosity. In short, the 
opinion that we derive from the previous discussion is 
that there is no serious problem with the data. Thus we 
move to interpret them within the 'standard scenario'. 

We analyzed the data of Kamiokande-II, 1MB and Bak- 
san adopting a definite model for neutrino emission . 
The first att emp t in this direction has been done by Lamb 
and Loredo [l7[. Of course, the choice of the model is a 
critical step. We assume the existence of the two main 
phases of antineutrino emission as expected in the 'stan- 
dard scenario'. For each phase we need to describe the 



1 This still leaves some open questions: The absolute times of 
Kamiokande-II and Baksan are not reliable. Is it fair to use 
the corresponding freedom of interpretation to set the first event 
in each detectors at t = 0? Do the assumed number of events 
fits into a reasonable model for the emission? Both of them are 
answered affirmatively with the global fit of the data discussed 
later, based on Il9l| . 




anti-v e energy [MeV] 

FIG. 5: The continuous curve is the V e spectrum in the best fit 
model; the dotted curve is the approximant described in Eq. [3 

luminosity of the emission, the average antineutrino en- 
ergy and the duration. For cooling, the parameters are 
R Cl T c and r c (similar to the parameters of the previous 
section). For accretion, we use the mass of neutrons 
exposed to the positrons, M a , the temperature of the 
positrons, T a , and the duration of the accretion, r a . In 
fact, the v e are produced from the reaction of the ther- 
mal positrons with the target neutrons around the proto 
neutron star, through e + n — > pD e , the inverse of IBD; 
a sample energy spectrum is shown in Fig. 2J Note the 
characteristic 'pinching' of the distribution-an output, 
not an input of our model. 

Let us mention for completeness some qualifying fea- 
ture of our analysis. We prescribe that (1) the tempera- 
ture of e + increases during accretion so that the average 
energy of v e is approximatively continuous at t ~ r a , 
that overcomes the shortcoming of the Lamb and Loredo 
model noted in [T3|; (2) the number of neutrons exposed 
to the positron flux decreases in time more smoothly than 
as in [171 ]: in this way the luminosity is also continuous, 
as expected on general basis. (3) Finally we avoid the si- 
multaneous presence of cooling and accretion v e delaying 
the cooling phase by an amount ~ r a again improving 
on [l7j. We also improved the energy spectrum of ac- 
cretion neutrinos and included neutrino oscillations with 
normal mass hierarchy. 

Our best fit result for the parameters is [191 ]: 

R c = 16 km, T c = 4.6 MeV, r c = 4.7 s, , . 
M a = 0.22 M Q , T a = 2.4 MeV, t q = 0.6 s. ( > 

The radiated binding energy is 8b = 2.2 x 10 53 erg, similar 
to the expected value for the formation of a neutron star. 
The values of the parameters of the cooling phase are 
reasonable; in particular, R c resembles a typical neutron 
star radius. Also the values of the parameters of accretion 
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R(0) 


P 


Ntot (z) 


Nthr (z) 


model 


2.00 
0.67 
1.25 


2.00 
3.44 
3.44 


3 0.5 
2 0.7 

4 0.7 


0.6 0.2 
0.3 0.2 
0.5 0.2 


ref.[23] 
ref.[13] 
ref.[24] 



TABLE I: 1 st and 2 nd column, parameters of SN distribution- 
R(0) is in units of 1/(10 4 Mpc 3 yr); 3 rd and 4 th column, total num- 
ber of expected events per year from the relic supernovae and aver- 
age redshift; 5 th and 6 th column, events above a threshold of 19.3 
Me V and average redshift. 



resemble the expectations; the value of initial accreting 
mass M a is a fraction of the outer core mass ~ 0.6M©. 
A last interesting outcome is that in the 30 s window 
Baksan had 0, 1, 2 background events with a posteriori 
probabilities of 20%, 47%, 29%, that compares well with 
the a priori expectation of 1 event. 

The electron antineutrino spectrum, shown in Fig. [5] 
can be calculated from the fluence (=time integrated 
flux) assuming D = 50 kpc, and it is well approximated 
by a modified Fermi-Dirac distribution: 



dN 
dE 



E- 



T 3 1 



exp 



' E 
. T 



with 



9.45 x 10 56 , 
= 3.77 MeV, 
0.531. 



(3) 

The a posteriori comparison with the detected energies 
is satisfactory, as already shown in Fig. [TJ 

As a possible application, we use this result to predict 
the number of relic neutrino events. We consider a de- 
tector a la SK with a fiducial mass of 22.5 kton (2^]. The 
rate of accumulation of IBD events is: 



N, 



dEa{E) 



Ho 



dz 



R( Z ) ^((i 



z)E) 



+ (W(l + zf 



(4) 



The cosmological parameters are Hq = 70 km/sec/Mpc, 
= 0.7, Qdm = 0.3. An important quantity is the rate 
of core collapse supernovae R as a function of the redshift 
z, adequately described by: R(z) = i?(0)(l + z) 13 if z < 1 
and = i?(0)2 /3 if z > 1. R{0) is the product of the fraction 
of core collapse supernovae fsN & n d the cosmic rate of 
star formation i?*(0). Other quantities (the maximum 
energy E v ~ 40 — 60 MeV and redshift ~ 5 — 6, the point 
z = 1 where the slope is modified) are not critical for 
the prediction. We compare in Tab.Uthe events per year 
for three descriptions of the rate of cosmic supernovae. 
The lower number of events is rather close to the one 
given in [l3j |: in other words, the non-standard neutrino 
spectrum used in [l3[ gives essentially the same results 
as the (more standard) spectrum of Eq. [3] Barring the 
possibility that SN1987A was a peculiar event, we are 
lead to believe that the largest uncertainty in predicting 
the relic neutrino signal is the rate of cosmic supernovae 
rather than the model for neutrino emission. 



IV. DISCUSSION 



The difficulty to interpret SN1987A data is due to 
small-number statistics, partial operativity of 1MB, poor 
absolute times in Kamiokande-II and Baksan, an ap- 
parent excess of events in Baksan, "gap" in the time- 
distribution of Kamiokande-II, peculiar angular distribu- 
tions and different energies of Kamiokande-II and 1MB. 
However, none of these problems constitutes an unsur- 
mountablc difficulty and on the contrary they can and 
have to be taken care in the analysis; furthermore, none 
of them can be solved too directly by considering rea- 
sonable deviations from the standard expectations. This 
does not guarantee us the possibility of reaching a safe 
interpretation of the data, but contributes to make less 
arbitrary the belief that these unique data can be under- 
stood in a coherent framework. 

A general question that we addressed is if the data re- 
semble the conventional expectations. We motivated the 
opinion that the presence of background events should 
be accounted for and that the signal should be described 
by considering an initial phase of intense neutrino lumi- 
nosity, as required in the 'standard scenario' HQ. This 
opinion is in agreement with the conclusions reached in 
1989 in the review on SN1987A by Imshennik and Nady- 
ozhin (H. 

Being aware of the incompleteness of the present in- 
formation (on 3D effects, rotation, magnetic fields; on 
peculiarities of SN1987A; on the existence of the neutron 
star; on neutrino oscillations; on the detailed detector 
response, etc [l|) we attempted to define the expecta- 
tions of the 'standard scenario' [1, 0] in a parameterized 
model, trying to keep it as simple and flexible as pos- 
sible but imposing certain general requirements, such as 
the continuity of the antineutrino luminosity and average 
energy. 

We elaborated on the pioneer ap pro ach of Lamb and 
Loredo to SN1987A data analysis [17| emphasizing its 
defects and drawbacks but arguing that their model for 
antineutrino emission can be improved to resemble more 
closely the expectations [l9j]. The numerical analysis 
of the data yielded indications in (perhaps surprising) 
agreement with the expectations of the 'standard sce- 
nario' for neutrino emission. We estimated in this way 
the spectrum of v e from the SN1987A data and ap- 
plied this result to evaluate the number of relic neutrinos 
events. 
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